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Abstract 

Recently evidence has been found that the perturbative QCD potential agrees well with 
phenomenological potentials and lattice computations of the QCD potential. We review 
the present status of the perturbative QCD potential and theoretical backgrounds. We also 
report our recent analysis which shows analytically, on the basis of renormalon dominance 
picture, that the perturbative QCD potential quickly "converges" to a Coulomb-plus-linear 
form. The Coulomb-plus-linear potential can be computed systematically as we include more 
terms of the perturbative series; up to three-loop running (our current best knowledge), it 
shows a convergence towards lattice results, e.g. At one-loop running, the linear potential is 
or with a = (27rCVA))A^. 



* Invited talk given at the "International Conference on Color Confinement and Hadrons in Quantum Chromo- 
dynamics (Confinement 2003)", Riken, Tokyo, 21-24 July 2003. 
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Figure 1: (a)Scale dependences of E to t(r) at r = 2tq ~ 1 fm, in the pole-mass and MS-mass schemes. 
(m = "T-Ms( m Ms) = ^ GeV) (b) Perturbative series of E tot (r) at r = 2r , and at /x = 0.53 GeV determined 
from dEtot/d/j, = in (a) (MS-mass scheme). 



1 Introduction 

In this article, we review the perturbative QCD predictions of the QCD potential for a static 
quark-antiquark (QQ) pair, in the distance region relevant to the bottomonium and charmonium 
states, namely, 0.5 GeV -1 (0.1 fm) < r < 5 GeV _1 (l fm). 

Several years ago, the perturbative prediction of the QCD potential became much more accu- 
rate in this region. There were two important developments: (1) The complete 0(a|) corrections 
to the QCD potential have been computed [1] ; also the relation between the quark pole mass and 
the MS mass has been computed up to 0(a|) [2]. (2) A renormalon cancellation was discovered [3] 
in the total energy of a static QQ pair*, E tot (r) = 2m po \ c + Vqcd(j)- Consequently, convergence 
of the perturbative series of E tot (r) improves drastically, if it is expressed in terms of the quark 
MS mass instead of the pole mass. 

Let us demonstrate the improvement of accuracy of the perturbative prediction for E tot (r) 
up to 0(a 3 s ), in the case without light quark flavors (nj = 0). In Fig. 1(a), we fix r to be 
2ro ~ 1 fm as an extreme long-distance case^ and show the renormalization-scale (/x) dependence 
of E tot (r = 2r ). We see that E tot is much less scale dependent when we use the MS mass instead 
of the pole mass. Fig. 1(b) compares the convergence behaviors of the perturbative series of E tot 
for the same r and when /x is fixed to the value where E tot becomes least sensitive to variation of 
/x (minimal-sensitivity scale). If we use the pole mass, the series is divergent, whereas for the MS 
mass, the series is convergent."'" In fact, we observe qualitatively the same features at different r 
and for different number of light quark flavors rx/. Generally, at smaller r, E to t{r) becomes less 
/x-dependent and more convergent. 

The aim of this paper is to study properties of _E tot (r), given the much more accurate prediction 

*By "static", we mean that the kinetic energies of Q and Q are neglected. 

Vq ~ 0.602 A^i [4] represents the Sommer scale. It is translated to physical unit following the convention of 
lattice calculations in the quenched approximation: ro ~ 0.5 fm. 

^When the series is convergent (MS-mass scheme), /x-dependence decreases as we include more terms of the 
perturbative series into E tot (r), whereas when the series is divergent (pole-mass scheme), xt-dependence does not 
decrease with increasing order. 
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Figure 2: (a) Comparison between Etot(r) and typical phenomenological potentials [6]. Arrows at the 
bottom show the r.m.s. radii of the heavy quarkonium states, (b) Comparison between E to t(r) and lattice 
computations of the QCD potential [7]. 



as compared to several years ago. In Sec. 2 we examine E tot (r) up to 0(a|). Sec. 3 provides our 
present theoretical understanding based on renormalon dominance picture and operator-product- 
expansion (OPE). We show that the perturbative QCD potential "converges" to a Coulomb+linear 
potential in Sec. 4. Conclusions are given in Sec. 5. 



2 E tot (r) up to 0(a|) 

The most solid part of our study is an examination of the perturbative predictions of E tot (r) up 
to 0(a|), through comparisons with phenomenological potentials and with lattice computations 
of the QCD potential. 

In Fig. 2(a), E tot (r) is compared with typical phenomenological potentials. Since the latter are 
determined from the heavy quarkonium spectra, we choose realistic values for the input parameters 
of E tot (r): mfe = 4.190 GeV; rii = 4 with m u = rn^ = m s = and fn c = 1.243 GeV for the light 
quarks in internal loops; fi = jj(r) is fixed by either dE tot /djj = or \E^\ = minimum, but both 
prescriptions lead to almost same values of E tot (r). See [5,6] for details. We see that E tot (r) 
corresponding to the present values of the strong coupling constant (dashed lines) agree well with 
the phenomenological potentials within estimated perturbative uncertainties (indicated by error 
bars). We also note that the agreement is lost quickly if we take as(Mz) outside of the present 
world-average values, so that the agreement is unlikely to be accidental. 

In Fig. 2(b), E tot (r) is compared with the recent lattice computations of the QCD potential 
in the quenched approximation. Accordingly we set n/ = in E tot (r). We take m = 2 GeV, 
which stabilizes the perturbative prediction up to largest r. The scale fi = fi(r) is fixed in the 
same way as in Fig. 2(a). See [7] for details. The three solid lines in Fig. 2(b) represent the same 
perturbative prediction for E tot (r)] there is an error in the relation between Aj^g and the lattice 
scale, and the three lines represent the error band of this relation Aj^gT = 0.602(48) as given 
by [4]. Taking this uncertainty into account and also uncertainties of the perturbative prediction 
(error bars), we observe a good agreement between E tot {r) and lattice results. 

In fact, by now several works have confirmed the agreements [5,8,6,9,7]. Although details 
depend on how the renormalon in the QCD potential is cancelled, qualitatively the same conclu- 
sions were drawn, i.e. perturbative predictions become accurate and agree with phenomenological 
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Figure 3: Schematic diagram showing n-dependence of |Vg"^(r)| (or l-^tcrb ( r ) I ^ n ^ ne pole-mass scheme) 
[black points] and that of l-E-tot ( r ) I i n the MS-mass scheme [red points]. 

potentials/lattice results up to much larger r than before. In particular, in the differences between 
the perturbative predictions and phenomenological potentials/lattice results, a linear potential of 
order AQ CD r at short distances r <C Aq^ d was ruled out numerically. 

3 Theoretical Backgrounds 

We would like to understand why the perturbative predictions of E tot (r) behave in the way we 
observed in the previous sections. In this section, we review our current theoretical understanding. 

3.1 Perturbative uncertainties: Leading-order renormalon 

The nature of the perturbative series of VqcdM and E tot (r), including their uncertainties, can 
be understood within the argument based on renormalons. This argument gives certain estimates 
of higher-order corrections in perturbative QCD, and empirically it gives good estimates even at 
relatively low orders of perturbative series. 

According to the renormalon argument, the perturbative series of Vqcd(^) (or E tot (r) in the 
pole-mass scheme) behaves as depicted in Fig. 3 [black points]; see e.g. [10]. Namely, the series 
shows apparent convergence up to order a 7 !? with n « 2tt/ (PqCcs)- Beyond this order, the series 
diverges. Due to the divergence (the series is an asymptotic series), there is a limitation to the 
achievable accuracy of the perturbative prediction. It can be estimated by the size of the terms 
around the minimum, 0(ag°); this gives an uncertainty for Vqcd(?") of order Aqcd ~ 300 MeV [11]. 

To obtain the above behavior, we may use, for instance, the "large-/?o approximation" [12] to 
estimate the higher-order corrections. The order ot^ +l term of VqcdO") for n ^> 1 is estimated as 

(B) ^ _2C £ aM 6/8 x rfasMy x n! (1) 

m ix V 2n / 

where (7p = 4/3 is a color factor and /3q = 11 — 2n//3 is the 1-loop coefficient of the beta function. 
From this expression, it is easy to obtain the behavior of the series as explained above. 

It is important to note that Eq. (1) is independent of r. Although Vp™\r) is dependent on r, 
the leading part in the large n limit is independent of r. This feature can be verified in Fig. 4(a): 
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Figure 4: V 0o {r) truncated at 0{a^ +1 ) term, En=o V ft W' We set ^ = 2 ' 49 GeV ' = 4 and " s ^) = °- 273 ' 
(a) Before cancellation of the leading-order renormalon. (b) After cancellation of the leading-order renormalon. 
Note that the vertical scales are different in two figures. (The figures are taken from [5].) 

As we include more terms of the perturbative series of Vqcd(^) (in the large- (3$ approximation), 
there are large corrections to the potential, which shift the potential almost r-independently. 
Furthermore, one sees no sign of convergence of the perturbative series in the range of r we are 
interested in. 

3.2 Renormalon cancellation and its implications 

As already stated, the leading-order renormalon uncertainty in the QCD potential gets cancelled 
in the QQ total energy, E tot {r) = 2m pole + Vqcd(?")> if we use the MS mass. As seen above, 
the leading-order renormalon is independent of r, hence, it can be absorbed by a redefinition of 
the quark mass. The pole mass is expressed in terms of the MS mass as a perturbative series, 
m P oic = m(l + Cicts + c 2 a| + ■■•). The leading part of this series at large orders takes the same 
form as Eq. (1) but in opposite sign and with half magnitude [13]. As a result, the leading-order 
renormalon is cancelled* in E tot (r) and its perturbative series behaves as [3] 



at large n. As compared to V/3 (r), the series converges faster and up to larger order, n « 
6tt/ (Po&s), but beyond this order again the series diverges; see Fig. 3. An uncertainty of this 
perturbative series can be estimated from the size of the terms around the minimum, which gives 
SE tot (r) ~ Aqcd x ( r AQ CD ) 2 . Thus, at distances r < Aq CD , the uncertainty becomes smaller 
than the original uncertainty <5Vqcd(?") ~ Aqcd- In fact, the error bars in Figs. 2(a,b) represent 
±|A 3 r 2 with A = 300 MeV. The agreement of E tot (r) and phenomenological potentials/lattice 
results holds within this uncertainty. 

A significant improvement of convergence can be seen by comparing Figs. 4(a) and (b), where 
higher-order corrections are estimated using the large-/?o approximation. Moreover, in Fig. 4(b) 
we see that the higher-order corrections make the potential steeper at large r, as compared to the 
tree-level (N = 0) Coulomb potential. In fact, due to this very feature, agreements between E tot (r) 
and phenomenological potentials/lattice results were observed in Sec. 2. We may understand the 

*For the complete cancellation of the leading behavior Eq. (f), one needs to expand Vqcd(?*) an d ?7i po i in 
the same coupling constant qs(/i). This is somewhat involved technically, since usually Vqct>(t) and m po i c arc 
expressed in terms of different coupling constants; see [5-7]. 




(2) 
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reason as follows. Let us define a strong coupling constant Of(//) from the interquark force: 



7-1/ \ d _ a F {l/r) 



(3) 



Since the leading-order renormalon in Vqcd(?") is ^-independent, it is killed upon differentiation by 
r. OLpil/r) grows at infrared (IR) due to the running (the first two coefficients of the /5-function 
are universal), which makes l-FfV)! stronger than the Coulomb force at large r. This means 
that Vqcd(^), after cancelling the leading-order renormalon, becomes steeper than the Coulomb 
potential. Within perturbative QCD, arguments based on F(r) are much more secure than those 
based on Vqcd(?") in the range of r of our interest. See [5,8] for details. 

Conventionally, theoretical calculations of the energy of a heavy QQ boundstate closely followed 
that of a QED boundstate such as positronium: it starts from the natural picture that, when an 
electron and a positron are at rest and far apart from each other, they tend to be free particles 
and the total energy of the system is given by the sum of the energies of the two particles (pole 
masses); as the electron and positron approach each other, the energy of the system decreases due 
to the negative potential energy, so that the total energy of the boundstate is given as the sum 
of the pole masses minus the binding energy. Applying the same description to the energy of a 
QQ system, however, is not natural, because when Q and Q are far apart from each other, the 
free quark picture is not good. As a result, the perturbative expansion of the boundstate energy 
turns out to be poorly convergent, due to the contributions from IR gluons with wave-lengths of 
order Aq^ d . On the other hand, intuitively we expect that there should be a way to calculate 
the boundstate energy in which the contributions of IR gluons can be mostly eliminated. This 
is because when the boundstate size is sufficiently smaller than Aq^ d , IR gluons cannot resolve 
the color charges of the constituent particles, so that they decouple from this color-singlet system. 
Such a calculation is possible if we use a quark mass (e.g. MS mass), into which only contributions 
from short wave-length gluons to the quark self-energy are absorbed (renormalized) . 

According to the renormalon argument (in the large- (3q approximation), the QCD potential 
and the pole mass are roughly given by 

VqcD(r) » - J e C F —^— , (4) 

m polc ~m + - J j^^f ~ 2 , (5) 

c/<m 

where q — \q\ and as(q) is (essentially) the one-loop running coupling constant. ^ It follows that 
qualitatively the QQ total energy can be expressed as [14, 5]* 

It shows that the energy is mainly composed of (i) the MS masses of Q and Q, and (ii) the self- 
energies of Q and Q originating from gluons whose wavelengths are shorter than the size of the 

'I'The large-/?o approximation is essentially a resummation of gluon vacuum polarization, which effectively changes 
the coupling constant ces(f-i) to the one- loop running coupling constant, as(q) = as(/u)/[l— (3oas(n) log(p,/q)/(2n)], 
where fi — e 5 / 6 /i in order to account for the non-logarithmic term of the vacuum polarization. 

■^IR part of Vqcd(0 in Eq. (4) is cancelled against that of 2m po i e , since e l $' r ~ 1 at g < 1/r. On the other 
hand, at g> 1/r, e l ^' r is highly oscillatory and contributions to Vqcd(^) from this region are suppressed. 
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Figure 5: Total energy of a heavy QQ pair is mainly composed of (i) the MS masses of Q and Q, and (ii) the 
contributions from gluons with wavelengths 1/m < X g < r. In (ii) the self-energies of Q and Q dominate over the 
potential energy between the two particles. 

system and longer than those absorbed into (i), i.e. 1/m < \ g < r. See Fig. 5. Thus, IR gluons 
(\ g > r) have decoupled from Eq. (6). If as(q) were a constant (non- running), E tot (r) would be a 
Coulomb potential. Due to the running of however, there is a large positive contribution in 

Eq. (6) as r increases and gets closer to A Q ^ D . This gives a microscopic description for how E tot (r) 
becomes steeper than the Coulomb potential at large r: It is the rapid growth of the self-energies 
(ii) as more IR gluons couple to this system with increasing r. — Again it stems from the running 
of the strong coupling constant. 

3.3 Consistency with operator-product-expansion 

An OPE of the QCD potential can be performed within the effective field theory "potential- 
NRQCD" [15]. The dynamical degrees of freedom in this effective theory have only very soft 
scale, where typical momenta are k <C 1/r. Higher momentum scales are integrated out and 
the effective interactions are generally non-local in space but local in time (potentials). In this 
framework, IR part of perturbative QCD calculations is factorized and absorbed into matrix 
elements of operators. A factorization scale /if ac is set to satisfy Aqcd <S /Ufac Cl/r. Then the 
total energy of a static QQ pair is given by [16,9] 

E tot (r) = 2m po i c + V PT (r; // fac ) + 5E vs (r; /x fac ), (7) 
5E VS = r 2 x^j^dt e - it( yo- y s) (^(t)0(t, 0) ab E b (0)) (/i fac ) + • • • (8) 

The perturbative potential Vpy(r; //f ac ) is now free from renormalons. The renormalons are ab- 
sorbed into 2m po i e and 5Eus(r; yUf ac ). In the leading-order of multipole expansion, 5E\js is pro- 
portional to r 2 , as shown above. This is consistent with our observation in Sec. 2. Namely, 
the agreement between E tot (r) and phenomenological potentials/lattice results holds within the 
0(Aq CD t 2 ) renormalon uncertainty, while a difference of order AQ CD r between them was ruled 
out (numerically) at distances r Aq^ d . According to OPE, indeed such a (linear) difference is 
not allowed theoretically, since there exists no operator that can absorb a linear potential. 

4 Pert. QCD Potential as a Coulomb+Linear Potential 

The observations and theoretical arguments in Sees. 2,3 indicate that the perturbative prediction 
of E tot (r) takes a Coulomb+linear form (up to an r-independent constant) in the range r < Aq^ d , 
as represented by a typical phenomenological potential. We would like to substantiate this claim 
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analytically. Since the uncertainty of the perturbative prediction is of order r 2 , there is a possibility 
that E tot (r) is predictable up to order r (linear) term at r < Aq^ d within perturbative QCD. 

Nevertheless, one immediately comes up with counter arguments: (1) A finite-order pertur- 
bative prediction of E tot (r) cannot be expressed in a Coulomb+linear form. This is because the 
perturbative QCD potential up to 0(ag ) has a form 

V N (r) = -C F ^^ xJ2®s(v) n Pn(log(ir), (9) 

n=0 

where P n {L) denotes an n-th degree polynomial of L. Therefore, as r — > oo, Vjv(r) — > 0, i.e. the 
tangent is zero, and there is no linear potential. (2) From dimensional analysis, the coefficient of a 
linear potential should be non-analytic in as, i.e. of order Aq CD ~ /i 2 exp[— 4ir/ ((3 a s )]- Therefore, 
it should vanish at any order of perturbative expansion. We will come back to these arguments 
later. 

Since we observed a numerical agreement between E tot (r) up to 0(a s ) and a Coulomb-plus- 
linear potential, it would be interesting to examine how E tot (r) up to 0(ag) or Vjv(r) will behave 
at large orders, N — > oo. For this analysis, we need (a) some estimates of higher-order terms of 
Vjv(r) and (b) some reasonable scale- fixing prescription. The latter is needed, since perturbative 
QCD does not provide any scale-fixing procedure by itself. Driven by such interests, recently 
we have shown analytically [17] that, based on the renormalon dominance picture, V/v(r) quickly 
"converges" to a Coulomb+linear potential for TV 1, up to an 0(Aq CD r 2 ) uncertainty. This 
was shown (a) by estimating the higher-order terms of Vjv(r) using renormalization group (RG), 
and (b) by adopting a scale-fixing prescription based on the renormalon dominance picture. 

Let us explain the scale- fixing prescription (b). According to the renormalon dominance pic- 
ture, if we choose a scale /i such that as(i^) ~ Qtt j ((3 N) , around this scale, Vn{t) (after cancelling 
the leading-order renomalon) becomes least /x-dependent and the perturbative series becomes most 
convergent; cf. Fig. 3. In view of this property, we fix \i such that* 

= K- 1 )- ( 10 ) 

Then we consider Vjv(r) for iV 3> 1 while keeping Aj^g finite. 

Let us first demonstrate our result in the simplest case. We set £ = 1 and the higher-order 
terms of Vjv(r) are estimated using the 1-loop running coupling constant [leading log (LL) ap- 
proximation] : 

Vn(t) = " / (0^ 4 ^ Mq)]„, (11) 

where [aizXoOK denotes the perturbative expansion of the 1-loop running coupling constant trun- 
cated at 0(ag): 



A? 



a s (n) 



MS 



«M E{^iog(^)}" (12) 



TV 



n=0 



Since the scale \x and as{p) are fixed at each iV via the relation (10), V/v(r) is a function of only 
r and once Aj^g is fixed. 

*Here, we generalize the prescription of [17] by introducing the parameter £, cf. [18]. 
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Figure 6: (a) v c (p), Cp and D(p, N) (N = 10, 30, 100) vs. p. (b) Comparison of v(p, N) — B(N) (N = 10, 30, 100) 
[dashed black] and the Coulomb+linear potential vc(p) + Cp [solid red]. (The latter is hardly distinguishable from 
the N = 30 line.) 



In the limit N ^> 1, Vjv(r) can be decomposed into four parts corresponding to {r l , r°, r 1 , 
r 2 } terms (with log corrections in the r _1 and r 2 terms): 

V N (r) = —^Av(Ar,N), (13) 
Po 

v(p, JV) = v c (p) + B(JV) + Cp + £>(p, N) + (terms that vanish as N -> oo), (14) 

where we rescaled r and V^v by A = p exp[— 27r/(/3 ai 1 g(p;))]. The "Coulomb"^" [t>c(p)], linear [Cp] 
and quadratic [-D(p, N)] parts are displayed in Fig. 6(a); see [17] for the formulas. The constant 
[B(N)] and quadratic [D(p,N)} parts are divergent as N — > oo, while the "Coulomb" [fc(p)] an d 
linear [Cp] parts are finite in this limit. The constant part diverges rapidly, but since it can be 
absorbed by a redefinition of quark mass in the total energy I£ tot (r), we will not be concerned 
with it. The quadratic part is divergent slowly as D ~ (p 2 /12) log N, and its size is small for 
r < A -1 and iV < 100 as compared to the Coulomb+linear part; see Fig. 6(a). As a result, when 
N is increased up to e.g. 10-30, v(p,N) quickly "converges" to vc(p) + Cp at r < A -1 , while it 
slowly varies as (p 2 /12) logiV, see Fig. 6(b). Note that, although v(p,N) in this figure have the 
form of Eq. (9), they approximate well t>c(p) + Cp at r < A -1 . 

If we vary £ in the scale-fixing prescription Eq. (10), as long as £ > 2/3, we obtain the 
same Coulomb+linear potential, vc{p) + Cp, whereas B(N) and D(p,N) change: for instance, if 
2/3 < £ < 1, D(p, N) is finite as iV — > oo, D{p) = p p x (log-corr.) with 1 < P < 2. 

We may consider that the constant [B(N)\ and quadratic [D(p, N)] parts represent renormalons 
in Vqcd(^)- If we perform OPE, these will be (partially) absorbed into 2m po \ e and 5E\js in 
Eq. (7). On the other hand, the Coulomb+linear part [vc(p) +Cp] may be regarded as a genuine 
perturbative part, which cannot be absorbed into 2m po \ e and 5E\js- 

The above result can be systematically improved by incorporating higher-loop effects [next-to- 
leading log (NLL), next-to- next-to-leading log (NNLL), • ■ ■] into the RG estimate of the higher- 

^The "Coulomb" part t>c(p) contains log corrections at short-distances and its short-distance behavior is con- 
sistent with the 1-loop RG equation for the potential. 
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Figure 7: Comparison of the Coulomb+linear potential (n/ = 0) with the lattice results in the quenched approx- 
imation. 1-, 2- and 3-loop correspond, respectively, to LL, NLL and NNLL approximations in the higher-order 
estimates. 



order terms of Vjv(r). Fig. 7 shows the Coulomb+linear potentials corresponding to the higher- 
order estimates by the 1-loop running (LL), 2-loop running with 1-loop non-logarithmic term 
(NLL), and 3-loop running with 2-loop non-logarithmic term (NNLL). They are compared with 
the lattice results. Since the 3-loop non-logarithmic term is not yet known, the 3-loop running 
case represents our current best knowledge. Up to this order, the Coulomb+linear potential agrees 
with the lattice results up to larger r as we increase the order. * 

The coefficient of the linear potential ar can be expressed analytically. For the 1- and 2-loop 
running cases, the expressions read 



2irC 

Cl-loop " 



AS° P ) 2 , (15) 



A, v MS J 

-oop - ^ (AS-) 2 r£r)[ 1 + a i s - 1 ~ Sei i {1 + * • (16) 

where 5 = (3\J (3$; see [17] for details. 

Finally, let us comment on the counter arguments given at the beginning of this section. (1) 
Viv(r) "converges" to a Coulomb+linear form for N 1. (See also discussion in [19].) In fact, 
already at relatively low orders and at r < A -1 , V/v( r ) approximates the Coulomb+linear potential 
fairly well. (2) Consider a perturbative term T = { \og(fir)} N . If we substitute the relation 
(10) and take the limit N — > oo, it is easy to see that T — - > (Ar) 3 ^. Thus, if we choose a scale 
/i a la renormalon dominance picture, perturbative terms can converge towards (Ar) p with some 
positive power P. Non-analyticity in as enters through the relation (10). 



5 Conclusions 

After cancellation of the leading-order renormalons by using the MS quark mass, perturbative 
uncertainty of E tot {r) = 2m po \ c + Vqcd(^) reduces from O(Aqcd) to 0(AQ CD r 2 ) and much more 

tThe NNLL originating from the ultra-soft scale [20] hardly changes the 3-loop running case displayed in 
Fig. 7 [17]. 
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accurate perturbative prediction is obtained at r < Aq CD ~ 1 fm. Consequently we observe the 
following: 

• E tot (r) up to agrees well with phenomenological potentials/lattice results within the 
estimated uncertainty. 

• New physical picture on the composition of the energy of a static QQ system or on the 
interquark force is obtained: The self-energies of Q and Q from gluons with \ g < r increase 
rapidly as r increases, which makes E tot (r) steeper than the Coulomb potential. 

• Our observations are consistent with OPE within potential-NRQCD framework. 

• Based on the renormalon dominance picture, the perturbative QCD potential "converges" 
quickly to a Coulomb+linear potential, which can be computed systematically as we include 
more terms in the estimate of higher-order terms via RG. In particular, the linear potential 
can be computed analytically [Eqs. (15,16)]. The Coulomb+linear potential agrees with 
lattice results up to larger r as we include more terms. 

The analyses reported in this paper can be applied to the bottomonium and charmonium 
spectroscopy in the frame of potential-NRQCD formalism, in which the QCD potential is identified 
with the leading potential in 1/m expansion. Some applications have already been done in our 
recent works [21]. Furthermore, our analyses provide a basis for (and justification for the error 
estimates of) the determination of the bottom and charm quark MS masses from the T(IS') and 
J/ip energy levels computed in perturbative QCD [14,22]. 
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